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ABSTRACT 


Based on the background of Carbon Peak and Carbon Neutralization (CPCN), this study aims to applies Remote Sensing 
(RS) images of Shenzhen in 2008, 2013 and 2018, combined with RS and Geography Information System (GIS) technology 
to classify Land Use/Land Cover (LULC), calculate Net Primary Productivity (NPP), and then estimate the carbon storage 
of green space in Shenzhen. The results show that during the decade from 2008 to 2018, the green space in Shenzhen is 
reduced and the construction land has increased. In the process of land transfer in Shenzhen, both cultivated land is non- 
agricultural and the reclamation of construction land is withdrawn, green space carbon reserves decreased in the 
beginning and then increased. The reason is that green space was transformed into other types of land from 2008 to 2013, 
resulting in the reduction of green space carbon storage. However, the change of green space from 2013 to 2018 is not 
obvious. Due to Shenzhen advocates a low-carbon economy and green development, resulting in an increase in carbon 


storage. 
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INTRODUCTION 


In the process of industrialization, greenhouse gas emissions dominated by carbon dioxide cause global warming and 
various accompanying extreme weather and natural disasters (Bonneuil et al., 2021; Cheng, et al., 2022; Fan and 
Wei,2022; Quevauviller,2022; Tao et al., 2022). Reducing greenhouse gas emissions has become the common 
responsibility and obligation of all countries. Facing a series of environmental problems caused by global climate change, 
the Paris Agreement clearly proposes to keep the global average temperature rise within 2 °C relative to the pre- 
industrialization level by the end of this century, and make efforts to control the global average temperature rise within 
1.5°C to reduce the risk and impact of climate change (Choudhury et al., 2022; Sun et al., 2022; Peng et al., 2022;Yu et al., 
2022). 


Following the trend of the times, China proposed in 2020 to strive for a national strategy to achieve the goal of 
carbon peak by 2030 and carbon neutralization by 2060(CPCN); (Qian, et al., 2022).Commonly, there are two ways to 
achieve carbon neutralization: one is to reduce carbon emissions, the other is to pay attention to the carbon sequestration of 
natural ecosystems (Li, 2021). Studies show that strengthening the carbon sequestration capacity of natural ecosystems is 


one of the most economical and effective ways to offset and absorb carbon emissions (Yang et al., 2022). 


With the development of China on urbanization policy, the control and monitoring of urban carbon emissions 
have a high correlation between CPCN (Shen et al., 2021; Imani et al., 2022; Tao et al., 2022). Due to the high degree of 
human activities and agglomeration in cities, besides causing a serials of urban heat island effect, urban air quality, urban 
ecosystem change, also they have a high carbon emission effect (Mahamba, et al., 2022). Therefore, how to pay attention 
to and solving the problem of urban carbon emission is an important indicator of government policy. Generally, green 
space, lakes, wetlands and other natural ecosystems have strong carbon sequestration capacity (Byomkesh et al., 2012; 
Dan-jumbo et al., 2018; Chen et al., 2019; Aka et al., 2022; Zhang et al., 2022). Besides contributing to the protection of 
urban ecology, its function also plays an important role in the cultural adjustment of living activities and leisure, and is also 
an important indicator of sustainable urban development. Thus, making good use of these carbon sequestration ecological 
spaces in the city will contribute to the monitoring and control of urban carbon emissions, and promote the policy of CPCN 


(here in after referred to as the “dual carbon” goal) (Qian, et al., 2022). 


Urban green space, in a broad sense, refers to all areas covered by vegetation within the city. As an important part 
of the urban ecosystem, urban green space can fix carbon and produce oxygen, and effectively alleviate the heat island 
effect. It is the main direct carbon sink way of the city. Carbon sink refers to the process, activity and mechanism of 
removing carbon dioxide from the air. Carbon sink in the terrestrial ecosystem mainly refers to the amount of carbon 
dioxide absorbed and stored by forests (Li, 2021), or the capacity of forests to absorb and store carbon dioxide. By 
enhancing the carbon sink function of urban green space, the city can achieve carbon and oxygen balance, or alleviate the 
carbon sink pressure of forest land outside the city, it helps to reduce the significant impact on production, life and ecology 
due to the increase of carbon dioxide concentration, and plays an important role in achieving the goal of "CPCN" (Shen et 
al., 2021; Yu and Wang, 2021;Qian, et al., 2022; Zhangand Zhou, 2021) Therefore, the quantitative study of carbon storage 


in urban green space has guiding theoretical and practical significance (Wang et al., 2022). 


As the first batch of low-carbon pilot cities, carbon emission trading pilot cities and sustainable development 
agenda Innovation Demonstration Zone (IDZ), Shenzhen has guiding significance for China's economic construction and 
development, low-carbon environment and green city. It is also the worthwhile study area of this paper. Thus, this study 
takes Shenzhen as the research area, uses Remote Sensing (RS) data and technology to realize LULC classification through 
Support Vector Machine (SVM) supervised classification method (Mirik et al., 2013; Chen, 2019; Talukdar et al., 2020; 
Xiong et al., 2020; Liou et al.,2021), meanwhile uses meteorological data as the data source to estimates Net Primary 
Productivity (NPP) through Carnegie-Ames-Stanford Approach (CASA) model(Potter et al., 1993 ; Field et al., 1995; Li et 


al., 2012), and finally using the vegetation-carbon-sequestration model to estimate the carbon storage of Shenzhen. 


This study is organized as follows: literature reviews are first presented, that refers to the main theory, concepts 
and technological trends in this study. In the following section, the study area and analysis methods as well as process 
framework are presents, along with data processing. The fourth and fifth section layout the analysis results of LULC 
change and the carbon store on the study area, then to figure out what’s real condition happen. Finally, a summary of the 


discovery on the monitoring benefits of this study and describing the contribution significance. 


LITERATURE REVIEW 
Carbon Storage of Urban Green Space 


As the most concentrated area of human life, cities account for a large part of carbon dioxide emissions. Meanwhile, with 
the acceleration of urbanization, carbon dioxide emissions are increasing. Studies show that 97% of carbon dioxide 
emissions come from urban areas. Thus, it is urgent to study the carbon sources and sinks of cities to provide decision- 
making basis for urban green and low-carbon development. Urban green space is one of the factors to increase urban 
carbon sequestration, so the quantitative research on the carbon sequestration capacity of urban green space is of great 
significance (Shen et al., 2021; Zhang et al., 2022). Generally, relevant research on carbon storage estimation of the forest 
ecosystem was carried out with the help of ecological methods (Guan et al, 1998; Kashaigili et al., 2013; Wang et al., 
2022).Carbon storage is the reserve amount of carbon element, the quality of carbon element or the amount of material, 
usually refers to the amount of carbon in a carbon pool (forest, ocean, land, etc.) (Wang et al., 2021; Wang et al., 2022).At 
first, the methods mainly include sample land inventory methods and biomass model methods. In terms of biomass also 
known as plant biomass, it refers to the dry weight of organic matter contained in unit area at a certain time. Later, it was 
gradually transformed into the measurement of carbon storage and carbon sink combined with RS technology (Yang et al., 
2005 ; Lu et al., 2008; Wang et al., 2021; Zhang et al., 2022). RS satellites and other technologies can dynamically acquire 
their change information in real time by means of space-time monitoring, which will play an efficient regulatory role in 


maintaining ecology and promoting carbon absorption. 


In the estimation of green space carbon reserves, the study is based on the estimation of forest carbon reserves 
first. After the measurement of forest biomass, the forest carbon reserves are converted at a certain proportion. The 
research on forest carbon reserves was first carried out in Germany, and then Japan, the United Kingdom and the United 
States gradually joined the ranks of carbon reserves estimation. Nowak (1993) estimated the total carbon reserves of 
American cities; Jo and McPherson (1995) estimated the carbon absorption and release of green space in Chicago 
residential area; Rowntree and Nowak (1991 )estimated the carbon storage of urban forests in the whole United States; 
Myeong et al. (2006) used RS technology to estimate the aboveground carbon reserves of urban forests in New York City 
by establishing a relationship model between NDVI and aboveground forest carbon reserves of urban forests. Bordoloi et 
al.(2022)used an integrated satellite-based approach to model space carbon stocks and carbon sink potential for different 


LULC patterns in northeast India. 


The research on carbon storage in China started to get involved in development in recent years. That is, with the 
continuous attention to global climate change, Chinese scholars have also started to study and estimate the carbon storage 
and carbon dioxide absorption benefits of urban forests. Such as Guan et al. (1998)outweigh the green space carbon storage 
in the built-up area of Guangzhou by using the measured data; Hanet al. (2003)introduced the application of RS and GIS 
technology in modern urban green space, and applied remote sensing and GIS technology to evaluate the ecological 
benefits of urban vegetation. Zhou et al. (2010) proposed an estimation model of carbon sequestration of urban green space 
vegetation based on remote sensing data in, and used the model to estimate the aboveground dry biomass and aboveground 


net primary productivity of urban green space vegetation. 


The NPP and CASA 


The Net Primary Productivity (NPP) of vegetation refers to the carbon absorbed by green plants through photosynthesis to 
remove carbon from their own respiration (Deng et al., 2022; Zhang et al., 2022). The remaining organic part can reflect 
the carbon fixation capacity of plants. It is of great significance for tracking regional carbon reserves, studying the carbon 
cycle, and helping to achieve “CPCN” (Xu et al., 2019; Zhang and Zhou, 2021). In terms of Carbon Peak, that is, the peak 
of carbon emission refers to the peak amount of carbon dioxide (or other major greenhouse gases) artificially discharged 
into the atmospheric environment at a certain time point (period). Additionally, Carbon Neutralization is also known as net 
zero carbon emission or net zero carbon footprint, its concept is divided into narrow and broad sense. Carbon neutralization 
refers to the net zero emission of carbon dioxide, but it is not the complete absence of carbon dioxide emissions, but the 
long-term balance between emissions and absorption. In a broad sense, carbon neutralization means not only carbon 


dioxide but also other major greenhouse gases (such as CHy, N2O, SF, etc.) should reach net zero emissions. 


As for the studies of NPP, the method of sample investigation and inventory was first used, focusing on the role of 
ecosystems in the carbon cycle with the method of biomass. With the development of time and technology, RS technology 
was introduced into relevant research, and RS technology can quickly and efficiently obtain the basic data of regional and 


large-scale ecosystems. 


In recent years, with the development of spatial analysis technologies such as RS and GIS (Amraoui et al., 2022; 
Uwiringiyimana and Choi, 2022; Zhao and Du, 2022), the light energy utilization model driven by RS data has developed 
rapidly, with CASA model as a typical representative. The combination of light energy utilization model represented by 
CASA model and RS technology provides convenience for the study of NPP model. The CASA model is widely used to 
estimate NPP at present. The CASA model builds an RS parameter model based on the mechanism of vegetation 
photosynthesis. After inputting meteorological data, land area, photosynthetic effective radiation of vegetation, coverage 
status and other parameters, it can estimate regional land NPP and global carbon cycle. It estimates NPP through RS data 
and meteorological data, which not only avoids the collection of complex parameters, but also can accurately simulate a 
wide range of NPP. In addition, the RS data adopted has strong periodicity, wide observation area and high spatial- 
temporal resolution, making it become a main method for estimating NPP and carbon reserves of forest ecosystem (Potter 


et al., 1993;Wang, 2021; Zhang et al., 2022) 


The research on the NPP of vegetation in China mainly began after 1970. Such as Piaoet al. (2003) estimated the 
NPP of vegetation in the Yangtze River Basin; Zhuet al. (2005) estimated the NPP of vegetation in Inner Mongolia, China; 
it is concluded that the carbon sequestration benefits produced by forests play an important role in the global carbon cycle; 
Liet al. (2015) estimated the carbon reserves of Shangri-La based on CASA model; Abdureyimu (2018) used the CASA 
model to inverse the NPP value of Ili Valley, calculated the vegetation carbon storage using empirical model, and analyzed 


the temporal and spatial distribution characteristics of vegetation carbon storage in Ili Valley. 
STUDY AREA AND METHODS 
Site Description 


Shenzhen, a prefecture level city under the jurisdiction of Guangdong Province, is a sub provincial city, specifically 


designated in the national plan, a megacity, and a special economic zone, a national economic center, an international city, 


a scientific and technological innovation center, a regional financial center, and a trade and logistics center approved by the 
State Council. Shenzhen has nine administrative districts and one new district: such as Luohu District, Futian District, 
Nanshan District, Baoan District, Longgang District, Yantian District, Longhua District, Pingshan District, Guangming 


district and Dapeng new district (show as Figure1). 


Shenzhen, a coastal city in southern China, is located in the south of the Tropic of cancer. It is located at 113° 46’ 
to 114° 37’E and 22° 27’ to 22° 52' N in the south of Guangdong Province, on the east bank of the Pearl River Estuary, 
adjacent to Daya Bay and Dapeng Bay in the East; It is adjacent to the Pearl River Estuary and Lingdingyang in the West. 
Shenzhen River is connected with Hong Kong in the South. It borders Dongguan and Huizhou in the north. Since the 40 
years of reform and opening up, Shenzhen has grown transform a small fishing village to an international metropolis step 
by step, from a traditional agricultural area become a global urbanization area, creating an urbanization miracle of 
"Shenzhen speed"(Wu et al., 2020).Becoming one of the four major central cities in Guangdong, Hong Kong, Macao, 
Dawan District, a national logistics hub, an international comprehensive transportation hub, an international science and 
technology industry innovation center, and one of China's three national financial centers, it will spare no effort be building 
a leading demonstration area of socialism with Chinese characteristics, a comprehensive national science center and a 


global marine center. 
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Figure 1: The Map of Shenzhen Administrative Division. 
Methods 


Determine the technical process of the research according to the main contents and methods of the research, hence this 
study proposes a processing framework (the technical framework shows Figure 2). The technical route is divided into four 
parts: such as data acquisition and preprocessing, classification of LULC types (including the land type transfer matrix 


calculation), NPP calculation as well as green space carbon storage calculation. 


e Data acquisition and preprocessing. The data in this study is using Landsat-8 satellite series data and 
meteorological data. The RS data are preprocessed by clipping, radiometric calibration, atmospheric correction 
and mosaic, and the solar radiation data are calculated by using the empirical model. Then, extracted the NDVI, 


Ratio vegetation index (SR), FPAR and APAR to preparing for the NPP of Vegetation calculation. 


e Classification of LULC types. Based on multi-source RS data fusion, the SVM classification method is used to 
classify the study area, and proving by kappa value, extract the regional LULC data, and then calculate the 


regional LULC dynamic degree as well as calculate the land type transfer matrix. 


e NPP calculation. Based on the classification results of urban LULC types and the results of physical and chemical 
parameters of vegetation, combined with meteorological data and the maximum solar energy utilization rate of 


vegetation, the CASA estimation model is established to realize the estimation of NPP. 


e Green space carbon storage calculation. Incorporation with NPP and LULC area data that were used to calculate 


the carbon storage of green space in Shenzhen. 
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Figure 2: The Framework of Process Methods. 
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Data Source and Processing 


The RS images used in this study are downloaded from the free imagery website which is provided by Geospatial Data 
Cloud (GDC). Landsat-8 images were downloaded in 2013 and 2018, but without landsat-8 images were available in 2008, 


so Landsat-5 images were used. The relevant data is shown in Table 1. 


Table 1: Metadata of Satellite Images used in this Research 


2008 | Landsat-5 TM 1,2,3,4,5,7 30m ino/dA 
2013 | Landsat-8 OLI 1,2,3,4,5,7,8 30m 122/044 
2018 | Landsat-8 OLI 1,2,3,4,5,7,8 30m 122/044 


This study uses ENVI 5.3 software to carry out a serial of process such as radiometric calibration, atmospheric 
correction, image mosaic and clipping for the images of the three periods respectively. The atmospheric corrected images 
are extracted, and the NDVI is calculated by near infrared (NIR) wavelength and red light band (R). Generally, the 
calculate expressions are used as follows: NDVI= NIR-R/NIR+R (Liou et al., 2021; Gouzile et al., 2022). 


Additionally, the NPP calculation involves solar radiation data. According to the existing empirical formula of 
total solar radiation (Du et al., 2003; Zhou et al., 2022), such as the solar radiation data of Shenzhen is calculated: Q = 
Q,(a + bs)WhereQ, is the total astronomical, radiationa, bare the empirical coefficients, ands is the percentage of 


sunshine. The sunshine hours are from the statistical yearbook of Shenzhen. 
Dynamic Degree of LULC 


LULC dynamics usually can describe the intensity of land change in a certain period of time and predict the trend of LULC 


change. The expressions are described as formula (1) (Li et al., 2021). 


Up—Va 


K= ay = =x 100% (1) 


Where: K is the dynamic degree of single LULC type in time; The interval is T; U, is the initial area of a land 
type study; U, is the area at the end of a land study. A positive value of K indicates an increase in land area and a negative 


value indicates a decrease in land area. 
Classification Method and Accuracy Evaluation 


Using ENVI software, fully considering the LULC characteristics of Shenzhen and referring to the land classification 
method of Maitiniyazi and Kasimu (2018), the LULC types of the study area are combined into five first-class LULC 
units, as shown in Table2, namely green space, water areas, cultivated land, construction land and unused land (Fatima and 


Javed, 2021), and then the SVM method used to supervise the classification to obtain the LULC data of Shenzhen. 


Table 2: The Classification of LULC pei Sl in Shenzhen 


imary Classification —__ 
Town 
Construction land Tend 
Settlements 
Industrial and mining land 
Reservoir 
Water areas Take 
Green space oo 
P Grassland 
Cultivated land Farmland 
Bare ground 
Unused land Oiier fands 


Visual interpretation is used to verify the accuracy of supervision and classification results in Shenzhen, and the 
overall accuracy is more than 90%. In this study, the overall classification accuracy and kappa coefficient are used to 
evaluate the accuracy of LULC classification results. The overall classification accuracy is equal to the sum of correctly 
classified pixels divided by the total number of pixels. The value of kappa coefficient indicates: 0.0-0.20 extremely low 
consistency, 0.21-0.40general consistency, 0.41-0.60 medium consistency, 0.61-0.80 high consistency and 0.81-1 almost 


completely consistent. 


NPP Calculation 


As an important part of surface carbon cycle, NPP is the main factor to determine ecosystem carbon sink. In the study of 
the impact of global change on the ecosystem, NPP has become an indispensable index (Tang et al., 2013).The principle of 
RS technology in carbon sequestration research is to calculate the relationship between absorbed radiation and primary 
productivity by recording the spectral response of plants during photosynthesis with sensors, so as to obtain NPP, and then 
estimate the carbon sequestration capacity of plants. According to the mechanism and structure of the model, the 
commonly used NPP estimation methods can be divided into three types: statistical model, light energy utilization model 
and process model. Based on the data of land classification results, in this study the improved CASA model is used to 
estimate the NPP of Shenzhen, gradually calculate NDVI, SR, FPAR and APAR, and finally calculate the NPP. The 


calculation process is as follows: 
(1) Normalized Vegetation Index (NDVI) 


In 1973, NDVI was proposed to observe the amount of vegetation on the earth's surface, which has become a widespread 
and practical method so far. The principle is that chlorophylls in plants has a large amount of absorption of red light (R) 
and strong reflection of near-infrared light (NIR). NDVI can be used to observe the vegetation cover of land. The closer 
NDVI is to 1, the higher the vegetation density is; the closer to 0, the lower the planting density. That is, the vegetation leaf 
surface has strong absorption characteristics in the visible red light band and strong reflection characteristics in the near- 
infrared band, which is the physical basis of quantitative RS monitoring of vegetation index (Liou et al., 2021; Gouzile et 


al., 2022; Zhou et al., 2022). 


Vegetation Index (VI) is a quantitative radiation value that reflects the relative abundance and activity of green 
living vegetation. It is often used to characterize the physiological status of vegetation, green biomass and vegetation 
productivity in the study area. At present, more than 20 kinds of VI have been proposed, such as Ratio Vegetation Index 
(RVI), Difference Vegetation Index (DVI), Normalized Difference Vegetation Index (NDVI), Enhanced Vegetation Index 
(EVI), Perpendicular Vegetation Index (PVI), Soil-Adjusted Vegetation Index (SAVI), Modified Soil Adjustment 
Vegetation Index (MSAVI), etc. Among them, the NDVI is the most widely used. The formula (2) for calculating this 


index is: 


NIR-R 
NIR+R ( ) 


NDVI = 


Where the NDVI range is [- 1, 1]. It is calculated by near infrared wavelength and red light band. 
(2) Ratio Vegetation Index (SR) 


The characteristics of plants reflecting and absorbing infrared light are used to reflect the richness of plants. However, this 
index will have large errors due to regional and seasonal differences, so it is impossible to directly use this index to 


compare the relative relationship of image planting in different regions or at different times. Represented by formula (3): 


1+NDVI (x,t) 


ce 1-NDVI (x2) 


Where NDVI (x, t) represents the NDVI value of pixel x in period t, which is calculated by the formula. 


(3) Fraction of Photosynthetically Active Radiation (FPAR) 


The FPAR refers to the effective radiation of photosynthesis incident above the vegetation canopy, that is, the absorption 
ratio of PAR (Photosynthetically Active Radiation) received by the vegetation canopy. Its calculation method is as follows: 
The maximum values (NDVI ,4x, SRmax) and minimum values (NDVI pin, SRmin) of NDVI and SR of vegetation types are 
calculated respectively through the cumulative frequency of 95% and 5%.Within a certain range, there is a linear 
relationship between FPAR and NDVI. FPAR can be determined according to the maximum and minimum values of 


NDVI of vegetation and the corresponding maximum and minimum values of FPAR(Zhu et al., 2006). 


The values of FPAR,,,,andFPAR,,;,are independent of vegetation type. 5% NDVI value represents desert and bare 
land, and 95% NDVI value represents that the vegetation is in full coverage. At this time, the corresponding FPAR, ax is 
0.95; 5% NDVI value represents bare ground, and its correspondingFPAR,,;, is 0.001(Zhu et al., 2006), which is expressed 
by formula (4): 


NDVI(x,t)—NDVIG min) 
NDVI max)—NDVI(i,min) 


FPARwpv1 = x (FPARmax — FPARmin) + FPARmin (4) 


The research shows that there is also a linear relationship between FPAR and ratio vegetation index (SR) (Zhu et 
al., 2006), which is expressed by (5): 


FPARcp = 2 Remi) (FPAR ae — FPARmin) + FPARmin (3) 
SR(imax)—SR(imin) 


In order to minimize the error, the study combines the FPAR estimated by NDVI and the FPAR estimated by SR, 
and takes the average value as the estimated value of FPAR. Since the non-vegetated area has NDVI value, the FPAR is 
calculated by region according to the current land classification situation of Shenzhen, ais the adjustment coefficient 
between the two methods. Take the average value of 0.5 (Maitiniyazi and Kasimu, 2018) and combine equation (4) and 
equation (5) to obtain equation (6): 


0 Construction land (6) 


EPAR(&,t) = { aFPARypy; + (1 —a@)FPARsp Green space 


(4) Absorbed Photosynthetically Active Radiation (APAR) 


The APAR refers to the proportion of photosynthetic effective radiation absorbed by vegetation in the total solar radiation, 
which reflects the strength of photosynthesis of vegetation, and the carbon sequestration capacity of vegetation. The light 
and effective radiation of vegetation depends on the total solar radiation and the characteristics of the plant itself (Jiao et 


al., 2020). Represented by equation (7): 
APAR(x,t) = SOL(x,t) x FPAR(x,t) x 0.5 (7) 
In formula (7), APAR(x, t)is the Photosynthetic active radiation (MJ.m™~.t') absorbed by the vegetation in 
pixel x in time period t; SOL(x, t)is the total solar radiation (MJ.m~.t') received by pixel x in time period t, which is 


calculated by empirical formula(Du et al., 2003); FPAR(x,t) is the proportion of photosynthetic effective radiation 


absorption of vegetation in pixels; 0.5 represents the proportion of total solar radiation used for photosynthesis. 


(5) NPP of Vegetation Calculation 
NPP of vegetation can be absorbed by plants, APAR and actual light utilization(e)expressed by two factors, the unit is (g 
C-m’): 
NPP(x,t) = APAR(x,t) x e(x,t) (8) 
€(x, t), represents the actual light energy utilization of pixel x in month t. 
Carbon Storage Calculation 


The following vegetation carbon sequestration model (Wang and Gu, 2012; Zhou et al., 2022) is mainly used to estimate 
the carbon storage of vegetation in the study area, which is estimated by the area of different land cover, NPP of vegetation 


and carbon conversion coefficient of different vegetation: 
E=Tx6xXc (9) 


Where, E represents the carbon sequestration amount of land cover type (kg); T refers to the land area (m°) 
corresponding to the land cover type; 6 represents the NPP of land cover type (gC-m’); c means the conversion factor 


between vegetation biomass and carbon content is 0.45 (Wang and Gu, 2012). 
ANALYSIS OF LULC CHANGE 


In this chapter, the spatial-temporal change analysis, dynamic degree analysis and transfer matrix analysis of LULC in the 
study area will be carried out respectively (Zhao and Du, 2022). Before this analysis process, RS images must be classified 
first, but the accuracy of classification is generally evaluated by kappa value. Finally, the kappa coefficient of the 


supervised classification results is above 0.90, indicated that with high accuracy. The results are shown in Table 3: 


Table 3: The Results of the Classification Accuracy Evaluation 


| 2008 | 2013 | 2018 
Accuracy Assessment 96.63% 96.46% 92.89% 
Kappa 0.95 0.95 0.90 


Analysis of Spatial-Temporal Change of LULC 


The supervised classification method selects the SVM to obtain the LULC data of Shenzhen in three periods, and then 
carries out classification statistics to obtain the spatial distribution area and area of various LULC types (as shown in 
Figure 3, Figure 4, Figure 5 and Table 4). In order to more intuitively unfold the LULC change in the study area, the 
LULC change area and LULC dynamic degree in Shenzhen were calculated. 
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Figure 3: The LULC Map of Shenzhen in 2008. 
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Figure 4: The LULC Map of Shenzhen in 2013. 
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Figure 5: The LULC Map of Shenzhen in 2018. 


Table 4: The Area of LULC Change in Shenzhen from 2008 to 2018 (unit: hm’) 
| 2008 2013 2018 


Land’ Type Area_| Proportion | Area __| Proportion | Area __| Proportion | 
Construction Land | 85490.7 43.8% 89056.0 45.6% 91886.2 47.0% 
Green Space 91677.6 46.9% 89531.3 45.8% 89178.4 45.6% 
Cultivated Land 2161 1.1% 3384.9 1.7% 2496.3 1.3% 
Water Areas 6997.5 3.6% 6925.3 3.5% 4854.2 2.5% 
Unused Land 9039.1 4.6% 6468.4 3.3% 6950.7 3.6% 
Total 195365.9 100.0% 195365.9 100.0% 195365.9 100.0% 


It can be found that the LULC types in Shenzhen are mainly construction land and green space, followed by water 
areas and unused land, and the cultivated land is the least (shown as Table 4). In 2018, for example, the construction land 
and green space area reached 91886.2 hm’ and 89178.4 hm”, accounting for 47.0% and 45.6% respectively, while the total 


proportion of cultivated land, water areas and unused land was only 7.6%. 


From 2008 to 2018, the construction land in Shenzhen increased by 6395.5hm*. As a carbon source, the 
construction of land is one of the important factors affecting carbon emissions. The vegetation is reduced, and the 
solidification effect of vegetation on carbon in the air is greatly weakened. Meanwhile, vegetation residues will also emit a 
large amount of carbon, and the solidification effect of soil on carbon will also be weakened, resulting in an increase in 
carbon emissions. From 2008 to 2013, the green area decreased steadily by 2146.3hm7, 2013 to 2018, which showed that 


Shenzhen attached importance to urban greening and low-carbon development. 


The cultivated land area increased slightly from 2008 to 2013, increased by 1223.9hm’, and decreased in 2018. 
The increase of urban population, rapid economic development and industrial structure adjustment are important reasons 
for the decrease of cultivated land. The water areas has decreased year by year, indicating that the development of 
urbanization has eroded a large number of water areas, which is not conducive to environmental protection in Shenzhen. 
The unused land showed a downward trend from 2008 to 2013 and increased in 2018, indicating that the land supply and 


demand in Shenzhen is tight, and the unused land is gradually developed for construction or other purposes. 


Dynamic Degree Analysis of LULC 


The single LULC dynamic degree mainly represents the change degree of a certain land type in a certain region. The 


dynamic degree of single LULC of Shenzhen in 2008-2013, 2013-2018 and 2008-2018 is calculated from equation (1). 


The results are shown in Table 5. 


Table 5: The Dynamics Degree of Single | LULC in Shenzhe 


Construction Land 


0.0083 % 


0.0064% 


0.0075 % 


Green Space -0.0047% -0.0008% -0.0027% 
Cultivated Land 0.1133% -0.0525% 0.0155% 
Water Areas -0.0021% -0.0598% -0.0306% 
Unused Land -0.0569% 0.0149% -0.0231% 


It can be seen from Table5 that during 2008-2018, cultivated land and construction land increased at the rate of 
0.0155% and 0.0075% respectively, while green space, water areas and unused land decreased at the rate of 0.0027%, 
0.0306% and 0.0231% respectively. The area of construction land increases year by year, with a rapid increase rate. The 
increase rate from 2008 to 2013 is faster than that from 2013 to 2018. The area of green space decreases year by year, but 
the reduction rate decreases from 0.0047% to 0.0008%. The single dynamic degree of the water areas decreased from 
0.0021% to 0.0598%. The reduction of the water areas will have an adverse impact on the ecological environment of 
Shenzhen. The rate of unused land has changed from decrease to increase, indicating that Shenzhen will develop other 


types of land to increase LULC. 
Analysis of LULC Transfer Matrix 


Applying the LULC area transfer matrix to quantitatively analyze the change direction of different land types and 
understand the evolution process of land types in the study area. It can be seen from Table 6 and Table 7 that the LULC 
change in Shenzhen during the two periods showed the characteristics of one land type changing to a variety of other types. 


Table 6: The LULC Type Transfer Matrix from 2008 to 2013 (unit: hm’) 


8696.8 1605.6 85490.7 


Construction Land 72484.4 1541.4 

Green Space 10207.9 76392.2 2014.8 2065.7 996.9 91677.5 
Cultivated Land 897.1 815.4 262.6 150.8 35.2 2161.1 
Water Areas 5906.6 1696.2 473.5 857.7 105.2 9039.2 
Unused Land 1351.2 341.7 81.4 322.4 4900.8 6997.5 
Total 90847.2 87942.3 3994.8 5002.2 7579.5 195365.9 


From 2008 to 2013, green space was the main type of transfer out, with a total of 15285.3hm’ transferred out, of which 
the area converted to construction land was the largest, accounting for 67.8% of the total transferred out area; Secondly, it is 
transformed into water areas and farmland. The decomposition of organic soil will be accelerated in the process of farmland 
cultivation, so the conversion from forest to farmland will lead to the release of greenhouse gases such as carbon dioxide from the 
terrestrial biosphere to the atmosphere. The construction land is mainly converted into green space, with 8696.8hm” transferred 
out. The withdrawal of construction land reclamation will increase the green area, which will increase regional carbon storage to 
a certain extent. The 897.1hm” of cultivated land is converted into construction land and 815.4hm? into green space. Cultivated 


land can be regarded as carbon sink to some extent. On the whole, good soil can store more carbon than release. 


Table 7: The LULC Type Transfer Matrix from 2013 to 2018 (unit: hm’) 


2018 
2013 Construction Green Cultivated Water Unused Total 
Land Space Land Areas Land 

Construction Land 80150.8 6032.2 460.9 1341.4 1070.7 89056.0 
Green Space 6218.8 80460.6 7719.6 1717.4 355.0 89531.3 
Cultivated Land 1135.4 1067.0 826.0 304.6 51.8 3384.9 
Water Areas 3459.0 1140.0 403.9 1323.0 142.4 6468.3 
Unused Land 922.3 478.7 25.9 167.8 5330.8 6925.4 
Total 91886.2 89178.4 2496.3 4854.2 6950.7 195365.9 


From 2013 to 2018, construction land and green space were still the main transfer out types. Due to population 
growth, transportation and housing will inevitably occupy a certain green space area. However, compared with 2008-2013, 
the transfer out area of green space decreased, indicating that Shenzhen attaches importance to protecting the ecological 
environment, and vigorously carries out to build the forest, to create a garden city, and strives to achieve the "National 
Forest City". Secondly, the conversion of cultivated land to construction land and green space is a positive ecological 
evolution process, which can significantly increase the carbon sink function of terrestrial ecosystem, which is also related 
to the implementation of the national policy of returning cultivated land to forest.Additionally, the water areasare mainly 


transformed into construction land. 


With the rapid development of Shenzhen, the transfer of construction land has gradually increased. Cultivated 
land, green space and unused land are the main sources of construction land. During this period, woodland has been 
transformed into cultivated land. It can be seen that before 2013, Shenzhen was mainly in a high-speed development 
period. In the process of urban construction, the demand for construction land is increasing, the scale of carbon source land 
is increasing, and the growth rate of construction land is large, it should lead to an increase in carbon emissions, and the 


forest land with high carbon sink is losing. 


NPP AND CARBON STORAGE ANALYSIS 
NPP Analysis 


According to the analysis and calculation of this study (as shown in Figure 6, Figure 7 and Figure 8), the maximum NPP of 
Shenzhen in 2008 was 160.6 g C-m’, the maximum value in 2013 was 134.0 g C-m’ and the maximum value in 2018 was 
134.4 g C-m’, showing the characteristics of high in the eastern mountainous areas and low in the central and western 
regions. Based on the analysis of the LULC classification map, the green space is mainly distributed in the economically 
backward southeast region, and the eastern region is mainly forest or mountainous region, with abundant vegetation types 
and relatively high vegetation NPP. The central and western regions have rapid economic development, frequent human 
activities, dense urban distribution, relatively small green space distribution and low vegetation coverage, so the NPP of 


vegetation is low. 


In the decade from 2008 to 2018, the NPP first decreased and then increased. From 2008 to 2013, Shenzhen was 
in the stage of rapid development. Obviously, construction land occupied a large amount of green space and destroyed the 
surface vegetation coverage, resulting in the decline of NPP. From 2013 to 2018, urbanization accelerated, a large amount 


of cultivated land was occupied, and the cultivated land area decreased. However, as can be seen from Figure 7 and Figure 


8, the NPP value of construction land in 2018 increased significantly compared with that in 2013. The reason is, during this 
period, Shenzhen also strengthened ecological construction, and increased the construction of road green space and park 
green space, resulting in an increase in the green area of Shenzhen. Thus, in 2018 Shenzhen won the title of "National 


Forest City". However, due to urban expansion, part of the green space was converted into construction land, so the NPP 


value changed slightly. 
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Figure 6: The NPP Distribution Map of Shenzhen in 2008. 
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Figure 7: The NPP Distribution Map of Shenzhen in 2013. 
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Figure 8: The NPP Distribution Map of Shenzhen in 2018. 


Carbon Storage Analysis 


According to formula (9), the calculation of carbon storage of green field in Shenzhen in 2008, 2013 and 2018 was 
4.26x10° t, 3.58x10° t, 4.2x10° t respectively. The results are shown in Table 8. 


Table 8: The Carbon Storage of Shenzhen in the Period 2008, 2013 and 2018 


Green area (hm*) 91677.6 89531.3 89178.4 
Greenland average NPP (g C-m’) 102.88 88.8 106.99 
Greenfield carbon storage (t) 4.26x10° 3.58x10° 4.2x10° 


The results show that the change of carbon storage of green space in Shenzhen first decreased and then increased. 
From 2008 to 2013, the change of LULC in Shenzhen led to the decrease of carbon storage of urban green space. During 
this period, the carbon sequestration capacity of urban ecosystem vegetation weakened. The main reason is that the green 
space with the highest carbon sequestration capacity was transformed into other types of land, and the green space area was 
transferred out by 15285.3 hm’. The decrease of the green space area led to a decrease in carbon storage of vegetation. 
However, from 2013 to 2018 the carbon reserves of urban green space are increasing, it is due to the fact that after 2012 the 
state has vigorously promoted the construction of ecological civilization and focused on promoting green development, 
sustainability development and low-carbon development. Driven by policies, the urban development in central and western 


regions of Shenzhen has gradually realized the transformation to green development. 


With the continuous expansion of urban construction in Shenzhen, the contradiction between land-use is 
prominent, and the high density of population and buildings makes the ecological environment face severe challenges. 


Therefore, Shenzhen has promoted a policy of "three-dimensional greening", which peculiar ideas has "borrow land" from 


roofs, walls, three-dimensional bridges, etc., and build greening in the air. This policy may also be one of the important 
factors affecting the change of carbon reserves of urban green space in Shenzhen. In this context, the ecological benefits of 


Shenzhen have been raised to a new level. 
CONCLUSIONS 


Based on the LULC change and NPP calculation of vegetation, this study estimates the carbon storage in Shenzhen and 


analyzes the change of carbon storage. The results show that: 


e During the period 2008-2013, the regional LULC type changes have obvious differences. Construction land 
increased by 3565.3hm”, cultivated land increased by 1223.9hm*’, and green space decreased by 2146.3hm. 
Among them, the decreased green space is mainly transferred out as construction land, and cultivated landas well 
as construction land are mainly transferred in. Thus, that kind of change we realize which resulting in the decrease 
of vegetation may also lead to the decrease of vegetation carbon storage, indicating that the carbon storage 


capacity of vegetation in Shenzhen was generally weakened during this period. 


e From 2013 to 2018, the difference of LULC types is not obvious, the change range of green space area is small, 
however the area of construction land is still increasing, so thatthe green space, and cultivated land, water areasas 
well as unused land are converted to construction land. Nevertheless, during this period, the NPP of urban green 
space showed an upward trend, meanwhile carbon reserves is increased by 0.95 compared with 2013 of 106 tons, 


which may has a certain relationship with returning farmland to forest and green policy development. 


From the perspective of above analysis results, this study adopts the combination of RS and GIS technology to 
analyze the land transfer of land environmental changes in Shenzhen in recent decade, as well as calculating the NPP of 
green space and carbon storage, which can effectively obtain the environmental monitoring benefits of the study area. 
Meanwhile, it can obtain further empirical proof according to the comparative evaluation of Shenzhen's environmental 
policies. Therefore, the experimental proof and contribution significance of this study have two crucial perspectives. 
Firstly, in terms of technical analysis, verify the feasibility and rationality of RS and GIS technology, as well as the 
availability and contribution in the long-term environmental perception monitoring in the future. Secondly, in terms of the 
practice of policy and decision-making, this study proves the effectiveness and contribution of urban greening, provides 
effective proof for the continuous promotion of the “dual carbon” policy in near future, and has a better empirical tool on 


the road of sustainable development. 


Therefore, to achieve the goal of “Carbon Peak and Carbon Neutralization”, we should vigorously develop urban 
greening, increase urban green vegetation, restore and protect forest, wetland, marine and other ecosystems, and increase 
the capacity of land and marine carbon sinks. In recent years, there have obvious implementing effectiveness on Shenzhen, 
that has actively built urban greenways, community greenways and regional greenways, and built a three-level park 
construction system of "Natural Park, Urban Park, Community Park", so as to form a green environment pattern with good 
ecology, beautiful landscape, distinctive characteristics and harmonious development between man and nature, and 


contribute to the realization of the environmental carbon neutral vision. 


Overall, the monitoring and analysis data of this study show that good greening policies and actions can indeed 


achieve the purpose of environmental improvement and provide feasible thinking for the realization of the “dual carbon” goal. 
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